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The ground state energies of Ag and Au in the face-centered cubic (FCC), body-centered cubic (BCC), simple 
cubic (SC) and the hypothetical diamond-like phase, and dimer were calculated as a function of bond length 
using density functional theory (DFT). These energies were then used to parameterize the many-body Gupta 
potential for Ag and Au. This parameterization over several phases of Ag and Au was performed to guarantee 
transferability of the potentials and to make them appropriate for studies of related nanostructures. Depending 
on the structure, the energetics of the surface atoms play a crucial role in determining the details of the nanos- 
tructure. The accuracy of the parameters was tested by performing a 2 ns MD simulation of a cluster of 55 Ag 
atoms - a well studied cluster of Ag, the most stable structure being the icosahedral one. Within this time scale, 
the initial FCC lattice was found to transform to the icosahedral structure at room temperature. The new set 
of parameters for Ag was then used in a temperature dependent atom-by-atom deposition of Ag nanoclusters 
of up to 1000 atoms. We find a deposition temperature of 500±50 K where low energy clusters are generated, 
suggesting an optimal annealing temperature of 500 K for Ag cluster synthesis. 



PACS numbers: 61.46.Df, 36.40.-c, 71.15.Mb,12.39.Pn 

I. INTRODUCTION 

Applications of noble metal nanoparticles are currently 
emerging in medicine where they are used as antimicrobial 
agents 1 , antivirals against HlV-fl, antiangiogenic agents^, in 
drug delivery 4 and in cancer therapy 45 . They are also used in 
chemical sensory devices due to their enhanced surface chem- 
ical activity 6 10 . Ag clusters as small as Ag7 and Ag§ sup- 
ported with mercaptosuccinic acid have been demonstrated 
to be very potent in water purification through their chemical 
sensitivity to the presence of heavy metals like Pb, Cd, Hg 11 . 
The application of noble metal nanoparticles in the electronic 
industry in inkjet printing of conductive lines for circuitory 12 
and as electronically conductive adhesivesS^EH results from 
their electronic properties. Organic memory devices based 
on DNA biopolymer nanocomposites with Ag nanoparticles 
have been demonstrated 17 . The large optical forces induced 
by the transfer of momentum from electromagnetic radiation 
to a dielectric nanoparticle also make them useful in nano op- 
tical manipulation l 7 l 18 lt 20 l. It has been shown that subjecting 
Au nanoparticles to optical forces induces aggregation and en- 
hances the Raman scattering intensity of the thiophenol cov- 
erage of the Au nanoparticle^. 

It is common to synthesize noble metal nanoparticle either 
as free standing or in an inert gas matrix with sizes rang- 
ing from one to several hundred of nanometers. Ag nan- 
oclusters as small as 4.1 nm and 5.6 nm have been deposited 
on silicon substrates. Their thermally induced disintegration 
has been studied, revealing melting at temperatures well be- 
low Ag melting temperature 2 ^. Clusters of few atoms are 
also routinely isolated by stabilizing them with some pro- 
tective molecules 22 2 3. For example, a 25 atom Au-thiolate 
cluster in solution 24 , dodecanethiol-stabilized-AuaJSl a 1 nm 
lyzozyme-stabilized-Au nanocluster for Hg 2+ sensor 2 ^ and 
DNA-encapsulated 10 atom Ag-cluster^have been reported. 
A small angle X-ray study on the mechanism of Ag nanopar- 



ticle formation showed that nanoparticle formation initiates 
with the formation of Agi3 clusters which agglomerate to- 
gether to form a nanoparticle within 6 mJ^l. 

The advances in synthesis methods have not been accom- 
panied by an equal increase in understanding of the struc- 
tural, electronic and optical properties. The nanosized nature 
gives these particles/clusters properties that are intermediate 
between molecular, which are of quantum origin, and bulk 
character. Size dependent structural changes have been widely 
investigated using methods ranging from experiments based 
on electron diffraction spectroscopy 2 ^! 3 -^ to theoretical meth- 
ods involving classical molecular mechanics approaches^ 
and quantum calculations. In particular, pseudopotential time- 
dependent density functional theory (TDDFT) calculations of 
the optical absorption of magic number noble metal nanopar- 
ticles have been carried oup3J Crossover from molecular to 
nanosized behavior was predicted to occur for clusters with 
more than 150 Ag atoms^l Calculations of the absorption and 
the Raman enhancement of Ag„ -pyridine complexes (n=2-20) 
have been done using TDDFT and both properties were found 
to depend strongly on the cluster size 3 -^. 

Thoretical studies of noble metal clusters are widely based 
on classical approaches making use of interatomic pote ntials 
like the Gupta potentiaP 12 ^, the Sutt on-Chen potentiaP 1221 
and the embedded atom potential 38 39 . The latter has been 
used to study the size dependent melting of Ag nanoclustet^S 
spontaneous alloying in Au-Ag nanoclusters 39 and structural 
optimization of Ag-Pd 41 . The Gupta potential is widely used 
in predicting stable structures of noble metal nanoparticles 
like Ag and Au 42 and bimetallic clusters like Pd-AvF 2 ^. The 
parameters of the potential are obtained by fitting experimen- 
tal or data on the bulk face-centered cubic (FCC) sys- 
tem. This does not, however, guarantee transferability of the 
potential to the different phases of the system. For example, 
when the resulting potentials are applied to lower coordination 
states of the material (such as dimers, trimers, etc.), very short 
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bonds and overbinding are obtained rendering the potential 
inappropriate for studies of low dimensional objects like nan- 
oclusters. This effect has also been demonstrated to explain 
the finding that the Gupta and Sutton-Chen parameterizations 
predict different growth patterns already for small Ag clus- 
ters^. Shao et al. 43 pointed out that due to the flatter nature 
of the dimer potential as given by the Gupta formulation, it is 
susceptible to yield more strain-tolerant structures than those 
generated using the Sutton-Chen potential. 

In this study, we use the DFT approach to perform ground 
state energy calculations on the FCC, BCC, simple cubic, the 
hypothetical diamond-like phase and dimer. These energy 
profiles are used to parameterize the Gupta potential. The pa- 
rameters are given as functions of the coordination number of 
the Ag and Au atoms-giving them a bond-order character. The 
appropriateness of the parameters for low coordinated struc- 
tures is demonstrated. 



II. ALL ELECTRON CALCULATION: 
PARAMETERIZATION OF THE GUPTA POTENTIAL 

We used the WIEN2k all-electron-full-potential linearized- 
augmented-plane-wave DFT code 4 ^ to calculate the binding 
energies of Ag and Au in FCC, BCC, simple cubic, diamond- 
like phases, and dimers. This DFT approach partitions the unit 
cell into muffin-tin (MT) spheres centered on each atomic site 
and the interstitial region. Within the MT spheres, the Kohn- 
Sham functions are given as linear combinations of spherical 
harmonics weighted by radial functions. In the interstitials 
they are parameterized as plane-waves. The generalized gra- 
dient approximation (GGApS for the exchange and correla- 
tion energy is adopted. The two parameters governing the ac- 
curacy of the calculation are the number of plane-waves in the 
Brillouin zone and the plane-wave vector k cut-off parameter 
- RKM. This latter parameter is the product of the maximum 
plane wave vector in the interstitial region and the smallest 
muffin-tin radius in the system. For our calculations, the total 
number of fc-points ranged from 1000 to 3000, depending on 
the system and RKM value of 6 were found to converge the 
binding energies to an accuracy of 3 mRy per atom. 

Non-magnetic calculations were performed since both bulk 
Ag and Au are known to be non-magnetic. For each of 
the noble metal structures, energy-bond length dependencies 
were obtained. From them, extrapolation to very large lat- 
tice parameters permitted the isolated atom values to be re- 
moved yielding the binding energies as shown in Figs [T] and 
[2] For FCC Ag, we find an equilibrium lattice parameter of 
4.07±0.02 A and a binding energy of 3.16±0.04 eV/atom. 
The former compares very well with the measured value of 
4.09 A while the binding energy overestimates the measured 
value of 2.95 eV^by about 7 %. The bulk modulus was found 
to be 106±2 GPa which is in excellent agreement with the ex- 
perimental value of 109 GPa at OK^S. For Au, a lattice pa- 
rameter of 4.09±0.02 A (see Fig. pla)) compares very well 
with the measured value of 4.08 AT The binding energy of 
4.01 ±0.04 eV, however, overestimates the experimental value 
of 3.81 eV 46 ' 4 ^' 49 '. The bulk modulus of 179±2 GPa is in ex- 



cellent agreement with experimentally measured value of 180 
GPaatOKP. 

The BCC phase is less stable than the FCC phase by about 
0.06 (0.08) eV/atom for Ag (Au) and is thus predicted to be 
the high pressure phase of these materials. This is in agree- 
ment with another theoretical work that predicted the same 
result for AvJSS. The small energy difference between the BCC 
and the FCC phases of Ag is a possible reason for the recently 
observed continuous and reversible BCC-FCC phase transfor- 
mation in Ag/V multilayer*^. 

We used these energy-lattice parameter (or energy-bond 
length E-r) dependencies to parameterize the Gupta potential 
for Ag and Au. For this, we rescaled the binding energies by 
a factor of 0.934 (0.951) for Ag (Au) to equate the FCC value 
with the measured value while rescaling distances by a factor 
of 1 .006 (0.996) for Ag (Au) to equate to the corresponding 
FCC experimental bond length of 2.889 A (2.884 A). With 
these, the bulk moduli become 99 GPa for Ag and 170 GPa 
for Au. Most of the parameterizations that have b een d one 
are based on the bulk properties in the FCC phas d 42 l 43 4 To 
our knowledge, this is the first attempt to develop parameters 
that make the potential transferable enough for use in various 
phases. These new parameters should therefore be appropri- 
ate in high temperature studies of these metals and for low- 
dimensional structures formed from them. 

For this potential, the energy of atom i is given by 52 



V, = v r - v a 

where the repulsive part of the potential is given by 



and the attractive part by 
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In the first approximation, the phase dependence of this po- 
tential is introduced by adopting a coordination dependence of 
the parameters; where coordination is defined by introducing 
a cut-off function f c 



, , v _ r 1 for x < 3.1 A 
U[X) ~\0 for x> 3.1 A 



by which coordination n ci is given by 

Wei — ^ fc{^ij/ 
j'eNei(i) 



(4) 



(5) 



where Nei(«) is the set of nearest neighbor atoms of atom i. 
The choice of a cut-off distance of 3. 1 A was made to exclude 
the next nearest neighbor atoms of the BCC phase that show 
up beyond 3.2 A for both Ag and Au. 
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FIG. 1. Energy per Ag atom as a function of distance between atoms 
for five phases of Ag: face-centered cubic, body-centerd cubic, sim- 
ple cubic, Ag in the diamond-like phase and Ag dimer. 
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FIG. 2. Energy per Au atom as a function of distance between atoms 
for five phases of Au: face-centered cubic, body-centerd cubic, sim- 
ple cubic, the diamond-like phase and Au dimer. 



The parameters A and /i are parameterized within the E-r 
dependence of FCC phases (Figs. [T]and [2]l as previous val- 
ues for Ag=2l and AvF^I did not permit a good fit of our DFT 
calculations. 

The parameters a(n c ), f3(n c ) and Ro(n c ) are given as func- 
tions of the coordination number of the atoms (and therefore 
of the phases). We suggest here the following simple func- 
tional forms: 



n c + n a 
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FIG. 3. DFT-parameters of the Gupta potential for the four phases of 
Ag: face-centered cubic, BCC, simple cubic, Ag in the diamond-like 
phase and Ag dimer. The fitted potentials are shown in panel (d) and 
the points are the corresponding DFT binding energies. 



where L(x) = 2 (exp (—a; 2 ) — 1 + a; 2 ) /x 4 . 

A fit of the parameters to these model functional forms are 
performed as shown in Figs.[3]and |4]for Ag and Au, respec- 
tively. In the fitting procedure, periodic images are included 
such that interactions between atoms separated by distances as 
far as 9 A are included. The parameters are listed in Table 

The effect of this new parameterization of Ag and Au clus- 
ters is that in lower coordination environments (such as sur- 
faces), the atoms become less energetic. For example, the en- 
ergy per atom of an Ag dimer was found to be 1 .09 e V. This 
compares better with the experimental value of 0.83±0.02 
eV ^1 tha n the previously obtained values of 1 .22 e V 4 ^ and 
1.4(PS] Also the dimer bond length of 2.59 A is in good 
agreement with the measured value of 2.5 3 APEED as com- 
pared to values of 2.37 A^and 2.43 A l 41 l 54 lthat result from the 
previous parameterizations. The best DFT dimer energy so far 
is 0.9 eV with a corresponding dimer bond length of 2.58 AE3. 
The Au dimer bond length of 2.56 A compares very well with 
the experimental value of 2.47 A^ whereas previous param- 
eterization sets this value at 2.3 AE2J. The dissociation energy 
per atom of 1.53 eV for Au2 is much closer to the experimen- 
tal value of 1.15 eV^than the value of 2.42 eV that previous 
parameterization gives 42 . The observations made here about 
the properties of dimers are also valid for other low coordina- 
tions such as 2-, 3- and 4-coordinated atoms. 

The implication of the weaker binding of lower coordinated 
atoms - as opposed to the strong binding obtained with the 
previous parameter^"^ - in studies of Ag and Au nanoclus- 
ters is that surface atoms will become more mobile. This re- 
sults in the nanoclusters experiencing surface premelting at 
lower temperatures. The current parameterizations may, thus, 
be very appropriate in surface studies. 

To check the validity of this parameterization for interac- 
tions beyond the cut-off distance of 3. 1 A as defined in Eq. [4] 
we considered the case where n c — and evaluated the pa- 
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TABLE I. Parameters of the Gupta potential for Ag and Au 
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FIG. 4. DFT-parameters of the Gupta potential for the five phases of 
Au: face-centered cubic, BCC, simple cubic, Au in the diamond-like 
phase and dimer. The fitted potentials are shown in panel (d) and the 
points are the corresponding DFT binding energies. 
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FIG. 5. The validity of the parameterization beyond r c =3.1 A. Single 
atom, two atoms and three atoms are clearly seen to approach the 
central Ag atom in such a way that the potential of the central atom 
closely resembles that of an atom in a dimer, a linear chain and in a 
trigonal arrangement of Ag atoms, respectively. 



rameters in Eq. |6| T he corresponding potential Uo(r) was ex- 
tracted from Eqs.|l]-[3] For this, the central atom is interacting 
with one, two or three equidistant atoms that are approaching 
it from beyond 3.1 A. In Fig. [5] we plot this potential for one, 
two and three atoms approaching the central atom. We also 
show the potential of the atom in the case where it participates 
in a dimer, in a linear chain and in a trigonal arrangement. The 
single atom approaches with a potential that is very similar to 
that of the dimer; two atoms approach with a potential very 
close to that of the atom in a linear chain while three atoms 
approach like trimers. The potential described here thus cor- 
rectly describes the short-, medium- and long-range interac- 
tions in Ag and Au. 



HI. RELAXATION OF AG55 CLUSTER AND MD 
DEPOSITION OF AG CLUSTERS 

We performed a 2 ns MD simulation to relax a 55 atom Ag 
cluster at 300 K. Temperature was kept constant by using the 
Nose-Hoover thermostat 6061 and the velocity-Verlet scheme 
was adopted for the MD moves. The cluster's structure was 
transformed from a 4-shell cluster to a 3-shell one within the 
first 10 ps. Figure [6] shows the density distribution as a func- 
tion of distance from the central atom. The cluster is charac- 
terized by a central atom, a Ag 12 first-shell of radius 2.74 A, a 
next shell of 30 Ag atoms of radius 4.73 A and an outermost 




FIG. 6. Left: Density distribution as a function of distance from the 
central atom of a Agss cluster for the unrelaxed (thin-dotted line) 
and the MD-relaxed structure at 300 K (full line). The feature at 
4.09 A in the unrelaxed cluster shifts upwards and merges with the 
feature at 5.0 A which shifts downwards. Right: Snapshot of MD- 
relaxed cluster at the end of the 2 ns simulation. 



shell of 12 atoms of radius 5.45 A. The structure of the cluster 
at the end of the 2 ns is shown in the right of Fig. [6] 

To test further our parameters, we performed an atom by 
atom molecular dynamic deposition of Ag nanoclusters of up 
to 1000 atoms. The starting configuration was an Ag 4 cluster 
with four Ag atoms sitting at the corners of a square. While 
the cluster is undergoing relaxation at a fixed temperature, an 
atom is initiated far away (4 nm from the center of the formed 
cluster) with a kinetic energy of 2 eV and velocity pointing 
toward the cluster center, with the position of the atom cho- 
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sen randomly on the sphere of radius 4 nm. Once this atom 
has entered the field of the relaxing cluster (distance less than 
3.5 A from the closest cluster atom) its dynamics are included 
in the Nose-Hoover thermostating scheme. The resulting clus- 
ter was relaxed for 100 ps and the process was repeated for 
subsequent atoms. A similar deposition process was done us- 
ing the old parameters 4 ^. 

Figure [7] shows the evolution of the energy per atom as the 
cluster size grows at different temperatures using both the new 
and old sets of parameters. We find that the energy decreases 
monotonically as the cluster size increases at all temperatures 
except at a temperature of 500±50 K, where sudden drops in 
energy were recorded for cluster sizes larger than 126 and 150 
for the old parameters, and 150 and 192 for the new set of 
parameters. The cluster size of 150 is equal to the reported 
size beyond which a cross-over from molecular to nanometer 
optical behavior has been reported^! We fitted the energy per 
atom as a function of the cluster size with the functional 

e(T) = E{N,T)/N = eoo (T) + bN~ 1/3 + cN- 2/3 , (7) 

where E(N, T) is the total energy of an N-atom Ag cluster at 
temperature T (in Kelvin). The temperature dependence of the 
extrapolated large cluster energy, eoo, is shown in the inset and 
reveals lower energy clusters at a temperature of 500±50 K. 
Note the big drop in energy as the temperature is raised from 
300 K to 500 K for the new parameters and the almost con- 
stant value of the energy when the old parameters are used. 
This indicates that the former parameters may be capable of 
distinguishing between the various temperature induced struc- 
tural changes within Ag nanoclusters. The finding of a more 
stable silver cluster at a temperature of 500±50K is in ac- 
cord with temperature programmed Auger, low energy elec- 
tron diffraction (LEED) and thermal desorption spectroscopy 
study of silver deposited on Mo(l 1 1) which demonstrated that 
stable 3 -dimensional silver clusters are formed once treated at 
temperatures above 300 K and below 650 KpJ. This tempera- 
ture range is also in agreement with the finding of an Ag crys- 
tallization transition between 250 and 310 °C where the ag- 
glomeration of silver atoms to form nanoclusters is favored 64 . 

Typical structural properties of a 561 atom Ag cluster atom- 
by-atom deposited at 500 K and annealed at 300 K with the 
new set of parameters is shown in Fig. [8] We calculated the 
electron diffraction pattern by using the Hartree-Fock atomic 
form factors as parameterized by Doyle and Turner 6 ^. The 
atomic distribution with respect to the central atom of the clus- 
ter G(r), shows a highly ordered nanocluster, with FCC order 
up to a radius of more than 1 nm. This is also confirmed by 
the structure factor S(q), with q = A being the wave 

length and 6 is the Bragg's diffraction angle. All results are 
averages over 10 ns in time intervals of 200 fs. The assign- 
ment of the spots agree excellently with that of Khan et al. 
on Ag nanoparticles 66 and Kang et al. on sintered inkjet- 
printed Ag nanoparticles^]. In particular, the electron diffrac- 
tion pattern is quite similar to the high resolution transmission 
electron microscopy (HRTEM) measurement of Ref. |66] on 
Ag nanoparticles of size 16.37 nm. G{r) also shows that the 
561 atoms cluster has a diameter of about 2.5 nm, which is 
also confirmed by size L of 24.97±0.21 A, obtained based on 
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FIG. 7. The evolution of the energy per atom vs Ag cluster size at 
different deposition temperatures. Open symbols: new parameters, 
closed symbols: olcF^. The arrows indicates points of energy drops. 
Inset: temperature dependence of (see Eqn.|7](. 



the position 2f?=38.16° and the full- width a t half-maximum 
(FWHM) cr=4.24° of the (1 1 1) spot a d 66 ' 68 ' 69 ' 



L 



0.94A 

a cos 28 ' 



(8) 



where a is in radians and the wavelength A=0.1541 nm is 
used. The 28 and a values are obtained from a fit (in the inter- 
val ±12° around the (111) spot) with a sum of three Gaussians 
and a Lorentzian, with the latter centered at q=Q. A similar 
calculation on a 2093 atoms cluster cut out of bulk Ag and 
relaxed at 300 K yielded a diameter of 4.22±0.29 nm which 
was also quite close to its size of 4.2 nm obtained on the basis 
ofG(r). 



IV. CONCLUSION 

We performed DFT calculations of the binding energies of 
Ag and Au in five different phases: FCC, BCC, simple cu- 
bic, diamond-like and dimer. We computed the corresponding 
equilibrium bond lengths, binding energies and bulk moduli, 
and found them to compare very well with experimentally 
reported values. We used the binding energy versus bond 
length dependencies to perform a cross-phase parameteriza- 
tion of the widely used many-body Gupta potential for Ag 
and Au. The new parameters, whose coordination-number 
dependencies are fitted to simple analytical functions, were 
found to correctly describe the energetic and structural behav- 
ior of low coordinated systems. We believe that these new 
parameterizations should be appropriate for studies of low di- 
mensional structures such as nanoclusters and surfaces. They 
may also be used to reconcile the structural differences re- 
ported for small clusters of noble metals generated using the 
traditional bulk-fcc-based parameters of the Gupta potential - 
which is susceptible to yield more strained structures - and 
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the Sutton-Chen potential, which leads to less-strained struc- 
tured. The current parameterization of the Gupta potential 
solves the problem of strain tolerance reported for the param- 
eters of Shao et alF^l as the broadness at the bottom of the 
dimer potential well is improved. 



40 60 80 100 120 
29 (decrees) 



FIG. 8. The lattice structure of the 561 atom Ag nanocluster gener- 
ated by MD deposition at 500 K and annealed at 300 K. The electron 
diffraction pattern (left) is calculated by using the Hartree-Fock pa- 
rameterized form factors of Doyle and Turner 65 . The (111) spot can 
clearly be seen. Other spots (200), (220), (311)/(222), (400) and 
(331)/420() are clearly visible and denoted as A, B, C, D and E, re- 
spectively. The density distribution with respect to the central atom, 
G(r), (top-right) shows a highly ordered crystal and the structure 
factor S(q) (lower-right) shows well defined features. The lines in 
the lower-right panel are for the functions used to fit the (111) spot. 
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